Deposition of Na Clusters on MgO(OOl) 
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We investigate the dynamics of deposition of small Na clusters on MgO(OOl) surface. A hier- 
archical modeling is used combining Quantum Mechanical with Molecular Mechanical (QM/MM) 
description. Full time-dependent density-functional theory is used for the cluster electrons while 
the substrate atoms are treated at a classical level. We consider Na6 and Nas at various impact 
energies. We analyze the dependence on cluster geometry, trends with impact energy, and energy 
balance. We compare the results with deposit on the much softer Ar(001) surface. 

PACS numbers: 34.50.Lf, 36.40.Sx, 68.49.Fg 
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I. INTRODUCTION 

A major branch of present-days cluster research com- 
prises clusters in contact with solid surfaces, for an 
overview see, e.g., (H-Q- The interaction of these 
two entities gives rise to a rich scenery of effects such 
as, e.g., chemical reactions at surfaces 0-0], particu- 
larly, catalytic applications (8l-Ho|. or modified optical 
response fllT - [l3l ] . One crucial aspect here is the process 
of cluster deposition which is relevant for synthesis and 
for analysis of clusters in contact with surfaces. More- 
over, the deposition dynamics as such is an interesting 
and demanding process due to the subtle interplay of 
the impact of interface energy, electronic band struc- 
ture of the substrate, and surface corrugation. Accord- 
ingly, there is a wealth of investigations on cluster deposi- 
tion, experimentally oriented |14T-(2ll]. theoretically with 
molecular dynamics (MD) techniques l22H26f or more de- 
tailed quantum mechanical methods [27H33I]. for reviews 
see @, H, |H[ . Although addressing the same physical 
processes, these various theoretical approaches, relying 
on different approximations, often provide useful com- 
plementary information. Recently, we have investi- 
gated deposition dynamics of Na clusters on Ar(001) sur- 
face [HI, IH, |36[ . The aim of this paper is to continue 
these theoretical studies now considering deposition dy- 
namics of Na clusters on a much "harder" surface than 
Ar, namely MgO(OOl) insulator surfaces. The struc- 
tural properties and optical response of Na„ on MgO 
have already been studied in great detail in using the 
present computational approach [33] ■ Both Ar and MgO 
are similar in that they are both insulators with a large 
band gap, but they differ significantly in other important 
properties. There are, however, large differences in other 
properties. Ar is a Van-der-Waals bound material, thus 
very soft with little surface corrugation. On the other 
hand, MgO is an ionic crystal, well bound and with large 
surface corrugation. It is thus most interesting to see 
how deposition dynamics proceeds in that such, 
and at variance with the Ar case. 

The theoretical description of clusters on surfaces is 



very involved due to the huge number of degrees-of- 
freedom of these systems. This holds the more so for 
dynamics. The vast majority of theoretical studies thus 
resorts to MD simulations using effective force fields be- 
tween the atoms, as mentioned above [2214261 . l3~ij . These 
are comparatively inexpensive and can provide a perti- 
nent picture of the leading atomic transport processes. 
Metal clusters are more than an ensemble of atoms be- 
cause they held together by delocalized bonds (due to a 
common electron cloud); in consequence they show pro- 
nounced shell effects [3814401 ] . This makes a quantum me- 
chanical description of cluster dynamics advisable. Most 
of the fully quantum mechanical pictures make compro- 
mises in concluding on dynamical features from a series 
of static calculations. A true Born-Oppenheimcr MD for 
deposition of Pd clusters on MgO substrate can be found 
in [29j]. The enormous expense of such high level cal- 
culations limits the size of the systems, particularly the 
size of the representative for the substrate. On the other 
hand, there are many situations in which the substrate 
is much more inert than the cluster. Our test case of Na 
clusters on MgO(OOl) belongs to that class. This sug- 
gests to use a hierarchical description where the cluster 
electrons are treated quantum-mechanically by full Time- 
Dependent Density-Functional theory (TDDFT) while 
the substrate atoms are handled at a lower level of refine- 
ment by classical motion. This modeling belongs to the 
family of coupled Quantum-Mechanical with Molecular- 
Mechanical methods (QM/MM) which are often used 
in other fields as, e.g., bio-chemistry [4l|-|43| or surface 
physics 0, [4f|. In earlier studies, we developed and ap- 
plied a QM/MM model for Na clusters in contact with 
Ar [H, Ha, 14614481 ] . We have shown that it was most cru- 
cial to include properly the dynamical polarizability of 
the substrate when exploring truly dynamical processes 
as we aim at. Recently, we extended the modeling to Na 
clusters on MgO surfaces, again including dynamical po- 
larizability [371 ] - Here we take up that model and apply 
it to a study of deposition dynamics. We will consider 
Na6 and Nas as test cases. These two clusters have very 
different geometries and binding properties which allows 
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TABLE I: The dynamical degrees-of-freedom of the model. 
Upper block : Na cluster. Lower block : Active cell of the 
MgO substrate. See text for details. 



to explore qualitatively the impact of cluster properties 
on the deposition process. 

The paper is organized as follows: In section [HI we 
summarize the QM/MM model for Na clusters on MgO. 
Section [Ml presents results tracking the detailed dynam- 
ics in terms of trajectories and analyzing the processes 
with respect to energy transfer and energy balance. Con- 
clusions are summarized in section HVl 



II. BRIEF SUMMARY OF THE MODEL 
A. The degrees-of-freedom 

The hierarchical QM/MM model has been detailed in 
[37l ]. We review that here briefly. The various con- 
stituents and their degrees-of-freedom are summarized 
in table III The Na cluster is treated in standard fash- 
ion [i^, [50| . Valence electrons are described in terms of 
single-particle wavefunctions ip n (r) and the complement- 
ing Na + ions are handled as charged classical point par- 
ticles characterized by their positions R^cn*) , see upper 
block of table D The electrons are described by TDDFT 
at the level of the local density approximation (LDA) The 
substrate is composed of two species: Mg 2+ cations and 
2 ~ anions. The cations are electrically inert and can 
be treated as charged point particles; they are labeled by 



The anions arc easily polarizablc, an aspect which is 



described by allowing for two constituents: a valence elec- 
tron distribution (labeled by i^) and the complementing 
core (labeled by i^). Each of these three types of con- 
stituents is described as a classical degree-of-freedom in 
terms of positions -Retype) , see the lower block of table |U 
The difference — represents the electrical dipole 
moment of the O 2- anion and is thus allowed, by con- 
struction, to explicitly evolve in time, as a function of the 
local electric field due to all constituents of the system 
at a given instant. 

The combined system is sorted in four stages of de- 
creasing activity, as sketched in figure [TJ The Na cluster 
is treated at the highest level of theory with full TDLDA- 
MD. The Mg and O ions of the substrate are arranged 
in fee crystalline order corresponding to bulk MgO, with 
a lattice parameter of 7.94 arj. All dynamical degrees- 
of-freedom for Mg and O, as listed in table HI are taken 
into account in an active cell of the MgO (001) surface re- 
gion underneath the Na cluster, denoted "zone I" in the 



sketch. The active cell is continued by an outer region of 
MgO material ( "zone Ha" ) where the ionic centers of Mg 
and O are kept fixed, while oxygen dipoles still remain ac- 
tive degrees-of-freedom. Thus zone I together with zone 
Ha constitute the "active cell". Anything farther out 
("zone lib") is totally frozen at crystalline configuration 
and only its Madelung potential is considered. The ef- 
fect of the outer region on the active part is given by 
a time-independent shell-model potential (5l| ; the actual 
parameters of this force field were adopted from EH] • The 
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FIG. 1: Schematic view of the hierarchical model for Najv on 
MgO(OOl) surface. 

active cell consists of three layers, each containing square 
arrangements of 242 Mg 2+ cations and 242 2 ~ anions. 
The ions in the lowest layer are fixed at the bulk structure 
to prevent them from relaxing and forming an artificial 
second surface. The volume where ions and electrons 
are mobile (zone I) has a diameter of 24 arj, all layers 
together (zone I+IIa) extend to 42 ao from the surface. 
Bulk structure farther out is modeled by the Madelung 
potential. Checks with models of a larger number of lay- 
ers showed that three layers provide an adequate descrip- 
tion of the very inert MgO material. (The soft Ar(001) 
substrate, used for comparison later on, is more critical 
and requires at least four active layers plus two frozen 
ones.) 



B. The energy 

The total energy is composed as E — i?Na + ^MgO + 
-Ecoupi where E^a describes an isolated Na cluster, i^MgO 
the MgO(OOl) substrate, and .Ecoupi the coupling be- 
tween the two subsystems. For E-^ & , we take the stan- 
dard TDLDA-MD functional as in previous studies of 
free clusters [49l. I50I ] including an average self-interaction 
correction [521 ] . The energy of the substrate and the cou- 
pling to the Na cluster consist of long-range Coulomb 
energy and some short-range repulsion which is modeled 
through effective local core-potentials [45j]. To avoid the 
Coulomb singularity and to simulate the finite extension 
of Mg 2+ and O 2- ions, we associate a smooth charge 
distribution p(f) oc exp (— r 2 /a 2 ) with each of these ionic 
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centers. We associate a similar smooth charge distribu- 
tion to the O 2- valence cloud as well. This altogether 
yields a soft Coulomb potential to be used for all active 
particles. 

C. Calibration of the QM/MM model 

The calibration of the whole model has to address three 
issues: The cluster as such, the environment as such, 
and the coupling between both. The modeling for the 
cluster is taken over from work on free clusters [H, H3] • 
The model parameters for the pure environment are the 
same as in previous studies of MgO(OOl) [45|, [Hj]. The 
parameters for the coupling between environment and 
Na cluster were calibrated from scratch. The tuning 
for Na@MgO(001) was performed using fully quantum- 
mechanically computed Born-Oppcnhcimer surfaces for 
Na atoms and Na + ions on MgO(OOl) from These 
surfaces were computed at four different substrate sites 
(0 2 ~, Mg 2+ , hollow, bridge) down to close distances 
where the full substrate repulsion was felt. For further 
details and the actual model parameters, see [13]. 

Two quantitative points are worth to be mentioned. 
The modeling achieves a barrier for penetration of clus- 
ter electrons into the substrate which reproduces nicely 
the large band gap of 6.9 eV for MgO. The fully quantum 
mechanical calculations show that electron transfer from 
the substrates O 2- anions to a Na atom remains below 
0.1 charge units down to the closest distances consid- 
ered (where the repulsive energy comes about the band 
gap). Transfer from the Na atom to the Mg 2+ cation is 
totally ignorable. This nice decoupling of ad-atoms and 
substrate is probably a feature of simple metals. Noble 
metals, e.g., can develop a more involved surface chem- 
istry due to the closeness of the d shell [2!| . 



pared to electronic motion. True electronic excitations 
are thus extremely small. For example, ionization stays 
safely below a fraction of 0.001 electrons. It would then 
be well justified to use Born-Oppenheimer-MD rather 
than full TDLDA-MD, as long as one carefully main- 
tains the crucial dipole polarizability of the substrate. 
But the TDLDA-MD scheme is so efficient that it is still 
preferable for reasons of computing time. Remind that 
the dipole response of the substrate needs to be propa- 
gated at electronic time scale and dipole stepping is more 
economic than fully relaxing the dipoles in each Born- 
Oppenheimer step. 



E. Preparation of the system 

First, the ground state structures of the pure MgO 
surface and of the free Na cluster are determined for the 
given model by simulated annealing. The Na cluster is 
then placed at a certain distance from the surface of the 
substrate. A distance of about 15 ao has turned out to be 
sufficient. After that a Galilean transformation is applied 
to the cluster. This means that each of the cluster ions 
is given a momentum Pq in the direction towards the 
surface and the electronic wave functions are boosted by 
an equivalent momentum as 

(pi(r) ->■ exp(zp • r) Pi(r) (1) 

where po = ^jj-Pq, iri e and M are the electron and Na 
ion mass, respectively. This provides the initial state 
from which on the system propagates in a straightforward 
manner according to the TDLDA-MD equations. 

F. Structure of the test cases 



D. Solution scheme 

From the energy functional, once established, one de- 
rives the static and dynamical equations variationally 
in a standard manner. The numerical solution of the 
coupled quantum-classical system proceeds as described 
in [37l I541 . l55l ]. The electronic wavefunctions and spa- 
tial fields are represented on a Cartesian grid in three- 
dimensional coordinate space. The numerical box em- 
ployed here has a size of (64 ao) 3 . The spatial derivatives 
are evaluated via fast Fourier transform. The ground 
state configurations were found by interlaced accelerated 
gradient iterations for the electronic wavefunctions [B(| 
and simulated annealing for the ions in the cluster and 
the substrate. Propagation is done by the time-splitting 
method for the electronic wavefunctions [5?J and by the 
velocity Vcrlct algorithm for the classical coordinates of 
Na + ions and MgO constituents. 

All the collisional processes studied in the current pa- 
per proceed on an ionic time scale, i.e. slow as com- 






FIG. 2: The structures of free Na 6 (left) and Na 8 (right). The 
bond distance along the five-fold ring of Nag is 6.5 ao and the 
top ion resides 3.1 ao above the ring. The bond distance in the 
two four-fold rings of Nas is 6.2 ao and the distance between 
the two rings is 5.8 ao. 

The starting point of deposition dynamics are well re- 
laxed structures for the clusters and pure MgO (001) sur- 
face. These had been discussed extensively in [37] . The 
MgO surface is a cut through cubic crystal structure. 
From the top, one sees a chess-board structure with al- 
ternating Mg and O ions. For the Na clusters, we will use 
here Na6 and Nas as examples. The initial state starts 
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from free clusters. Their structures are shown in figure [21 
Note that the vertical axis in the figure will represent the 
direction perpendicular to the surface in the forthcom- 
ing deposition processes [z axis). Nag is strongly oblate 
consisting out of a ring of five ions topped by one sin- 
gle ion. Nag has a highly symmetric configuration out of 
two rings of each four ions tilted relative to each other by 
45° to minimize Coulomb energy. The electronic cloud of 
Nag is close to spherical shape because N = 8 electrons 
correspond to a strong shell closure for Na clusters [loT ]. 
It is important to note that the bond distances for Nag 
(6.2 oo) are not far from the diagonal distance between 
oxygen sites in the MgO(OOl) surface (5.7 do) while the 
dimensions of the fivefold ring in Na6 do not fit well to 
the surface. That will play a role in the dynamical evo- 
lution studied later on. The equilibrium distance of the 
lower cluster plane (facing towards the surface) and the 
first surface layer is 5 oq. 



III. RESULTS AND DISCUSSION 

A. Na monomer on MgO — the influence of sites 

A two component system like MgO has more possible 
adsorption sites than a homogeneous material like an ar- 
gon substrate. The properties of an oxygen site are much 
different from those of the magnesium site because of the 
much larger polarizability of oxygen. We will thus con- 
sider four positions with respect to the surface : O site, 
M g si te, hollow and bridge. The structure calculations 
of [33] have shown that O, due to its large polarizability, 
is the most attractive site while Mg acts like a repulsive 
site on the cluster. This is in agreement with quantum 
chemical ground state calculations of transition metals 
on MgO [Hi] . 

In order to check the adsorption properties of the var- 
ious sites, we first study the deposition dynamics of a 
Na monomer. We briefly remind the static properties 
of Na@MgO(001). The O site is most attractive, bind- 
ing Na 5 a above the surface with energy 0.25 eV. The 
Mg site is dominantly repulsive. The hollow and bridge 
sites lie in between these extremes. For deposition dy- 
namics, the atom was initialized 15 ao away from the 
substrate above an O site, Mg site or hollow site respec- 
tively, each with an initial momentum along z direction, 
pointing perpendicular towards the surface with a mag- 
nitude corresponding to a kinetic energy E® ia = 0.136 
cV. Figure [3] shows the results of the simulation. The 
^-coordinates are chosen such that the (average) MgO 
surface layer resides at z = 0. The simplest case is the 
impact on the O site. The atom approaches the surface 
up to a distance of 4.5 a which is reached at about 500 fs 
and transfers part of its momentum to the substrate ions. 
The transfer proceeds at a very short time scale. The 
surface itself is excited mainly by the first collision which 
initially only affects the ions in the immediate vicinity 
of the atom at closest impact. The perturbation quickly 



spreads over the surface, but the associated sound wave 
does not penetrate very deep into the surface. The os- 
cillations in the third layer are already almost negligible. 
After the instant of closest contact, the atom bounces 
back, but it has already lost so much energy that it can- 
not escape from the surface anymore. Thus it performs 
damped oscillations, with each bounce transferring some 
momentum to the surface and being practically adsorbed 
within the first 2ps. The final distance approaches nicely 
the equilibrium distance of 5 ao. The right upper column 
of figure [3J shows the corresponding kinetic energy con- 
tributions. In the first 400 fs, the attraction from the 
MgO substrate leads to a rapid increase of the kinetic 
energy of the atom up to 0.45 eV. At the point of clos- 
est contact, the repulsive part of the interface potential 
stops the atom abruptly. That first collision transfers by 
far the largest amount of energy to MgO, whereas at all 
subsequent collisions the energy decreases more slowly. 
In order to check that the oscillations proceed only per- 
pendicular to the surface, the kinetic energy of the Na 
atom has been split into contributions from perpendicu- 
lar (or vertical) and parallel (or transverse) motion. The 
latter is too small to be visible in figure [3J and practi- 
cally negligible. Thus the motion of Na proceeds strictly 
perpendicular to the surface. The kinetic energy trans- 
ferred to the MgO can also be read off from figure [3J 
(see right upper panel). The contributions from oxygen 
and magnesium are given separately. Oxygen ions are 
the lighter species and therefore react first being quickly 
accelerated. About 100 fs later, the energy has already 
been distributed almost equally over both ion types. 

The dynamics behaves totally different if the atom im- 
pinges on the (repulsive) Mg site, see middle panels of 
figure El At first glance, the z-component of the Na tra- 
jectory looks quite similar to the case before. But one 
notes that the motion is not damped after the first re- 
flection. The kinetic energies (middle right panel) give a 
clue on the process. There is much less energy transfer 
at first impact which is related to the fact that the Mg 2+ 
ion is more inert. And there is a significant amount of 
lateral kinetic energy for the Na atom creeping up after 
impact time at 500 fs. In fact, most of the kinetic energy 
is now in lateral motion. The atom is deflected by the 
Mg 2+ ion. It is to be noted that the annealing of the 
substrate configuration leaves a small amount of symme- 
try breaking with fluctuations of the atomic positions of 
about 0.05 ao- This small symmetry breaking allows the 
atom to acquire sidewards momentum and so it bounces 
away in sideward direction, hops over the surface sev- 
eral times changing direction whenever it comes close to 
another surface ion. The motion is almost undamped be- 
cause little energy is transferred to the surface after the 
first collision. The atom has thus still too much energy 
to be caught by a certain site of the surface. But as the 
atom cannot escape the surface as a whole, it will con- 
tinue to lose slowly energy and finally be attached to an 
oxygen site, long after the simulation time of 3 ps. 

The bottom panels of figure [3] show the case of impact 
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FIG. 3: Time evolution of the ionic coordinates and kinetic energies Ekin for the deposition of a Na monomer on MgO(OOl) 
with initial kinetic energy E^ in = 0.136 eV, impinging on various sites : O site (top), Mg site (center) and hollow site (bottom). 
Left : z coordinates of Na (thick line), Mg (thin curve), and O (gray line) cores. Right : Total Eki n of Na (thick gray or red 
curve), lateral Ekin of Na (magenta dots), vertical -Ekin of Na (black dashes), and total Ekin of Mg (thin dark or blue curve) 
and O (thin light or green curve) cores. 



at an hollow site. We see again the immediate reflection 
at impact time associated with fast energy transfer. Less 
energy is transferred than on the other sites (see upper 
and middle panels) and thus the bounce-back has a much 
larger amplitude than in both other cases. The Na mo- 
tion remains strictly perpendicular to the surface as prac- 
tically no lateral kinetic energy can be seen. The vertical 
kinetic energy is almost approaching zero because the 
departing Na atom has to work against the polarization 
potential. The case is at the limits of our box size and 
energy resolution such that we cannot decide whether the 
atom will finally escape with extremely small kinetic en- 
ergy, or will bounce back and relax to an adsorption site 



on a very long time scale. Nevertheless, we find it worth- 
noting that the hollow site seems sufficiently attractive to 
hinder deflection towards the still more attractive oxygen 
site. 



B. Cluster deposition 

1. The case of symmetric Na% 

The analysis of section llH Al has shown the importance 
of surface site nature in the deposition process. Deposit- 
ing an extended object such as a cluster will lead to a 
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time (ps) ionic x-coord. (a ) 

FIG. 4: Left : Time evolution of the ionic ^-coordinates of 
Nas approaching MgO(OOl). Right: Projection of the Na and 
MgO trajectories into the x-y plane; the figures are vertically 
sorted by increasing initial kinetic energy : Soft deposition 
(top), robust deposition (center), reflection (bottom). 



mixed situation because the ions of the cluster will nec- 
essarily be placed above different sites. In the following, 
we will discuss deposition of Nas an d Na6 which have 
very different structures and so promise to show different 
deposition scenarios. As a first step in the analysis, we 
shall consider detailed ionic trajectories both perpendic- 
ular and parallel to the surface. The case of Na§ is shown 
in figure |4j The cluster was injected with its symmetry 
axis pointing through a hollow site and with the lower 
ring facing closer to bridge sites. The top panels show a 
soft deposition where the initial kinetic energy is 0.109 eV 
(0.0136 eV per Na ion). The left upper panel shows that 
the cluster is slightly accelerated in the initial phase, due 
to the attraction from the surface. But that acceleration 
differs for the different ions on the lower ring because 
they approach different sites on the surface. At the same 
time, the cluster rotates in the x-y plane to bring the 
four ions of the lower ring closer to the attractive oxy- 
gen sites. One may spot that from the top view in the 
right upper panel. At the point of closest impact around 
900 fs, the cluster transfers some momentum to the sur- 
face. The substrate ions are slightly displaced from their 
equilibrium positions and oscillate around their new po- 
sitions. The disturbance quickly decreases from layer to 
layer. The perturbation is negligible already in the fourth 
layer. The complicated detailed dynamics of the Na ions 
indicates that a major part of the translational kinetic 
energy is converted into heat, i.e. kinetic energy of the 
motion relative to the center of mass, as will be confirmed 



in section UlI Dl Nevertheless, the cluster basically keeps 
its original structure during the whole simulation period 
of 9 ps. In particular the two rings, each made of four 
ions, always stay clearly separated from each other. The 
top-down projection of the trajectories (see right-hand 
side of figure H|) shows that the cluster as a whole (or 
its center of mass) remains oscillating around the point 
of impact. The remaining kinetic energy of the cluster 
does apparently not suffice to overcome the surface cor- 
rugation barriers. This is related to the fact that the 
Na§ structure fits approximately well to the structure 
and binding distance of MgO, see section Hi Fl 

The middle panels of figure |4] show a more robust 
deposition dynamics with initial kinetic energy E®- ln = 
1.09 eV. The initial velocity is higher and the impact 
time comes earlier, now at 450 fs. The pattern remains, 
in principle, similar to the softer deposition. There is lit- 
tle momentum transfer to the substrate, strong internal 
excitation of the cluster, and the cluster is not departing 
too far from the impact point. However, perturbations 
are much larger, yielding larger amplitudes in vertical 
and lateral motion. As a consequence, the two rings of 
Nas are now not always clearly separated. Nevertheless 
the typical structure of Nas reappears from time to time 
as we will see later. The top-down projection of the tra- 
jectories (middle right panel of figure indicates a new 
effect, a sideward drift from one adsorption site to the 
next equivalent site. This sideward drift is again induced 
by the interplay between attractive O and repulsive Mg 
sites. The chaotically moving Na ions explore a strongly 
corrugated surface which leads to occasional side kicks 
from the repulsive Mg sites. 

The bottom panels of figure 2] show a hard collision 
with initial kinetic energy E^ in = 10.9 eV. Internal clus- 
ter and excitation and surface perturbation are, of course, 
again larger. The new feature is that the cluster is re- 
flected from the surface and leaves the numerical box at 
about 1 ps, however with huge internal excitation. It is 
not clear whether the departing cluster will stay asymp- 
totically stable. That is beyond our simulation capacity. 

Figure [5] complements the view within showing a se- 
quence of snapshots of the detailed structure for each of 
the three cases discussed above. The uppermost panel 
for soft deposition nicely shows the initial rotation of 
the cluster to match the attractive oxygen sites. The 
further snapshots indicate the sizeable internal excita- 
tion, however remaining small enough to see at all times 
clearly the two-ring structure of Na 8 . The middle panel 
for more robust deposition also presents the much larger 
cluster oscillations where the original cluster structure is 
often blurred, but reappears shortly at other times. That 
demonstrates the surprisingly good binding of Na clus- 
ters, particularly the Nag cluster with its magic electron 
configuration. The lowest panel shows the case of re- 
flection. Obviously, some ions would like to stick to the 
surface, but are finally caught back by the cluster which 
departs in a highly excited state. 
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FIG. 5: Snapshots of the time evolution for three cases of col- 
lisions of Nag on MgO(OOl) with different initial kinetic ener- 
gies i?kin- Upper panel : Soft deposition with E^ in = 0.109 eV. 
Middle panel : Robust deposition with E® in = 1.09 eV. Lower 
panel : Reflection with E° in — 10.9 eV. 



2. The case of strongly oblate Na® 

Results for the deposition of Na6 are shown in figure 
[5] The impact energy is varied and all three cases start 
from the same initial configuration where the top ion of 
Nag (see section III F|) is facing away from the substrate 
and the fivefold ring is parallel to the surface. The gen- 
eral features are similar to the case of Na§ . One observes 
a large internal excitation of the cluster while compara- 
tively little perturbation goes to the substrate and there 
is again the clear distinction between deposition for lower 
impact energies and reflection for higher ones. But there 
are several interesting differences in detail. Most of all, 
there is a strong lateral drift in all cases. Indeed the 
pentagonal ring of Nag does not match the rectangular 
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FIG. 6: Left : Time evolution of the ionic z-coordinates of 
Na6 approaching MgO(OOl). Right : Projection of the Na and 
MgO trajectories into the x-y-plane; the figures are vertically 
sorted by increasing initial kinetic energy : Soft deposition 
(top), robust deposition (center), reflection (top). 



structure of MgO, which hinders it from fully accomodat- 
ing the attractive oxygen sites. Thus one or two corners 
of the pentagon are bent up during the deposition pro- 
cess, and this, in turn, induces a sizeable lateral momen- 
tum (see right panels of figure |6]), and a strong perturba- 
tion of the pentagon, as can be deduced from the motion 
of z-coordinates shown in the left panels. In the case 
of the robust deposition, the cluster even rolls over the 
surface. The stronger lateral excitation leaves also some- 
what more perturbation to the substrate than in the case 
of Nas , as may be spotted when comparing figures [5] and 
^ That will become more obvious when checking ener- 
gies in section lHIDI Finally, it is interesting to note that 
the cluster orientation is also reverted in the case of re- 
flection. The z coordinates (upper left panel in figure [5]) 
suggest a process where the ring and the former top ion 
are reflected "independently" such that the topping ion 
is departing "behind/after" the ring, and thus reverting 
the cluster orientation. 

As noted above, when deposited, the Na 6 experiences 
a sizable drift due to the mismatch of its structure with 
the crystalline structure of MgO. One thus expects that 
direction and strength of the lateral motion depend sen- 
sitively on the initial position and orientation of Nag rel- 
ative to the surface. Indeed figure [7] shows the trajec- 
tory of the center-of-mass of the Nag cluster projected 
onto the x-y-plane for three different initial orientations. 
There are obvisouly dramatic differences. The cluster is 
kicked to a strong lateral motion for initial impact at 
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x[a ) 

FIG. 7: The trajectories of the Na6 center of mass, projected 
on a plane parallel to the surface. The lateral drift sensitively 
depends on the site above which the cluster impinges (hollow, 
O or Mg site) and on its orientation (top ion above - up - 
and top ion below - down - the ring). 



slow deposition process (left panels) , there is a strong in- 
crease of the z-polarization at the time of closest impact. 
This polarization remains during the further evolution 
at about 10% of the Wigner-Seitz radius, hence repre- 
sents a considerable internal polarization. Also, some 
x-y polarization builds up during the ongoing deposition 
oscillations. In contrast, in the case of a reflection (right 
panels), one sees a large instantaneous polarization at the 
time of closest approach, but only a very small remain- 
ing effect when the cluster has departed from the surface. 
The very short interaction time limits the internal exci- 
tation. 



C. MgO versus Ar Substrate 



repulsive sites (Mg, hollow) while only moderate lateral 
drift appears for impact on the attractive O site. 

One can learn more about the electronic charge distri- 
bution in the cluster during a collision by a direct multi- 
pole analysis of this distribution. We discuss here briefly 
the lowest non-trivial moments, the dipoles. Figure [8] 
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FIG. 8: Dipole polarization of the cluster electrons during 
collision of Na6 on MgO (001). Upper panels: Dipole moments 
of the electron cloud of the Nag cluster in x-y- (horizontal) 
and z-direction as function of time. Lower panels: Time- 
evolution of the z componente of the center-of-mass of the 
Na6 cluster; the line z = corresponds to the center of mass 
of a fully relaxed configuration of Na6 on MgO(OOl). Left 
panels: case of soft deposition at rather low initial kinetic 
energy (as indicated). Right panels: case of more energetic 
collision which leads to immediate reflection. 

shows the time evolution of dipole polarization for two 
different scenarios, deposition vs. reflection. The figure 
is augmented by the time-evolution of the center-of-mass 
as global indicator of the dynamical situation (lower pan- 
els). The two scenarios differ by the initial kinetic ener- 
gies, the lower value related to a more or less soft depo- 
sition, while the higher initial velocity leads to immedi- 
ate reflection of the cluster, as discussed above. In the 



In previous works |33j, |4]]], the deposition of Na6 on 
a cold, condensed argon substrate Ar(001) was investi- 
gated. Like MgO, solid Ar has a large band gap. But 
apart from this insulating nature, the two materials have 
much different properties. The attractive interaction be- 
tween MgO and Na is much stronger than between Ar 
and Na, due to the larger polarizability of the oxygen 
ion. On the other hand, frozen Ar material is a very 
soft solid due to the weak Ar-Ar binding. The melting 
point of Ar is 83.78 K (5t|, much lower than that of MgO 
around 3073 K • The softness of the Ar material thus 
changes the energy balance to the extent that the Ar 
substrate takes up most of the impinging energy, leaving 
rather little internal excitation for the cluster itself. Thus 
Ar substrates are very efficient soft stopper materials. In 
the former analysis [33|, S3] , we had run a similar series of 
impact energies as above and we also found soft deposi- 
tion for energies up to at least E® in /N- mn = 0.272 cV. An 
attempt to reach a reflection regime by further increasing 
the impact energy then led to a significant destruction of 
the substrate. Figured] illustrates that violent collision of 
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FIG. 9: Time evolution of z coordinates during collision of 
Na6 with initial total kinetic energy E^ n = 8.16 eV (1.36 eV 
per Na atom) on an Ar(001) surface. Data taken from [47j . 

Na6 on an Ar(001) surface. The cluster is indeed finally 
reflected, but the process evolves much different from the 
case of the collision with MgO shown in figure |U The 
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cluster is not reflected instantaneously, as for MgO, but 
with a delay of about 500 fs. It requires an additional 
boost from momentum reflected by the first Ar layer to 
finally release the cluster. Moreover, so much energy has 
been deposited in the weakly bound Ar material that 
the substrate is seriously damaged by that forced "reflec- 
tion" . These significant differences between Ar and MgO 
substrate will also be seen in the energy analysis later on. 



D. Energy Transfer 

1. Time evolution of energy components 
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FIG. 10: Time evolution of various contributions to kinetic 
energy for the collision of Nag with MgO(OOl) at three differ- 
ent initial kinetic energies as indicated. Contributions are : 
Total kinetic energy of the cluster ions (Na, thick curves), con- 
tributions due to center of mass motion (Na transl., dashes) 
and relative motion or heat (Na relat., dots), and total kinetic 
energy of the substrate (MgO, thin curves). 



into centcr-of-mass energy and intrinsic kinetic energy 
(from the motion relative to the center-of-mass). Let us 
first consider the case of reflection (upper panel). In the 
approaching phase, the cluster is accelerated by about 
0.68 eV which is small compared to the initial energy 
10.9 eV. Dramatic and fast changes emerge at impact 
time at 200 fs. The cluster kinetic energy exhibits a deep 
minimum. In that stage, almost all energy is stored in 
deformation. A large part of that deformation energy 
is quickly released showing up now as intrinsic kinetic 
energy of the cluster plus a smaller bit in translational 
energy. Another small fraction of energy is transferred 
to the substrate. The translational kinetic energy de- 
creases further on, because the departing cluster has to 
work against the attractive polarization interaction. Still, 
there remains sufficient translational energy to allow the 
cluster to finally escape, as in a very inelastic collision. 

Similar results are found for the soft (lower panel in 
figure [Til)) and robust deposition (middle panel). Again, 
only a small fraction of the energy is transferred to 
the substrate, another small fraction goes to the cluster 
ccnter-of-mass oscillations, and the major part to intrin- 
sic energy of the cluster. Particularly interesting is the 
case of robust deposition (middle panel) where it requires 
a second bounce to stir up intrinsic cluster motion. Be- 
fore that, there is still enough energy in translation to 
allow a lateral hopping from one attractive MgO site to 
the next (see also figured]). It is worthnoting that the 
average trend of the kinetic energy of MgO in figure [TU] 
has a small, but nonvanishing, slope. The cluster con- 
tinues to exchange energy with the substrate on a very 
slow pace. That indicates a thcrmalization process which 
eventually leads to equidistribution of kinetic energies af- 
ter long time, however much beyond our simulation ca- 
pabilities. 



Na 6 with MgO(001) 

E kin =8.16eV 
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FIG. 11: Time evolution of the dipole polarization energy 
in the MgO (001) substrate for the collision for Na6 with 
MgO(OOl) at an initial kinetic energy of 8.16 eV (reflection 
case). 



A complementing view of the deposition dynamics is 
given by the kinetic energies. Figure [10] shows the time 
evolution of the kinetic energies for Na and MgO. The 
kinetic energy for the Na cluster is furthermore splittcd 



The present modeling includes an independent dynam- 
ics of the dipole moments of the oxygen anions in the 
substrate. It has been shown recently that a signifi- 
cant amount of energy can be stored in these degrees- 
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of-freedom when a metal cluster is deposited on an Ar 
surface [IB, HfJ- One can define a dipolc energy which 
scales as the square of the dipole amplitudes. Figure [TT] 
shows a typical result for the time evolution of the en- 
ergy contained in the oscillating dipolcs in the case of 
reflection of Na6 deposited on MgO. There is a small 
initial value which corresponds well to the finite initial 
distance of the Na cluster to the substrate. There is a 
large contribution at the time of closest impact. This 
part is dominated by (instantaneous) static polarization 
which would also be contained in a Born-Oppenheimer 
MD. The dipole energy falls back to lower values when 
the cluster departs from the substrate (see figure EJ). But 
there remains some offset which corresponds to the en- 
ergy finally transferred to the dipole degrees of freedom. 
It amounts to about 2% of the impact energy, which is 
small compared to the other energetic obscrvables, see 
figure [14] and corresponding discussion in section IIII D 31 
In contrast to the case of Ar substrate, energy transfer 
to MgO dipoles has actually only a small effect for the 
overall ionic dynamics. But more subtle properties as 
optical response and trajectories of free charges (to be 
discussed in a subsequent publication) will be sensitive 
to such details. 



2. Energy transfers "at" impact 

Notwithstanding asymptotic thermalization, the fast 
energy transfer to the substrate in the early stages is 
an interesting observable characterizing the collision pro- 
cess. Figure [T2l shows the kinetic energy of the substrate 




Ekin (eV) 

FIG. 12: Heating of the MgO or Ar substrate after impact 
of Na6 and Nas. The energy transfer to Ar is independent of 
the cluster structure, whereas the energy transfer to MgO is 
not. Na6 transfers about twice as much energy as the compact 
Na 8 . 

soon after the collision, i.e. averaged over the first 2 ps 
after impact, as a function of the initial kinetic energy of 
the cluster E^ in . Apparently the energy absorbed by the 
substrate is proportional to E® in . But the slope depends 
very much on cluster and surface types. The soft Ar sub- 
strate absorbs much more energy than MgO, typically a 
bit more than 50% of the initial kinetic energy. The 
softness and the rather small surface corrugation of Ar 



make the process insensitive to the actual cluster which 
is approaching. That is different for MgO. There is al- 
ways less energy absorption by the substrate and there 
is a strong dependence on the cluster configuration. Nag 
transfers more than twice as much energy as Na§. The 
strong surface corrugation of MgO induces that sensitiv- 
ity to cluster geometry Remind that Nag does not match 
very well to the MgO surface while Nag does (see section 
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FIG. 13: Upper panel : The gain in kinetic energy of the clus- 
ter when approaching the MgO surface; lower panel : Instan- 
taneous energy loss, defined as the difference between max- 
imum kinetic energy before the impact and next maximum 
after impact. 

More information about what is happening directly 
around impact time can be obtained by reading off ob- 
scrvables at shorter time scales (shorter than the 2 ps 
used above) . Figure [T3] shows two such observables as a 
function of initial kinetic energy. The upper panel shows 
the energy gain in the approaching phase due to the ac- 
celeration by the polarization potential. It is defined as 
the difference of the first maximum of the cluster kinetic 
energy and the initial energy. Nag acquires slighlty more 
energy than Nas because it has a non-vanishing dipole 
moment which, in turn, enhances the polarization attrac- 
tion. There is a large gain in the low energy range (the 
regime of soft deposit), while the trend becomes very 
flat for fast collisions. This is probably due to a move 
from adiabatic to non-adiabatic relaxation processes in 
the surface. For very low impact velocities, the sur- 
face ions have time to follow the forces from the cluster, 
whereas for very high velocities the surfaces ions do not 
have enough time to respond before the cluster collides. 
The lower panel of figure [T3] shows an attempt to quantify 
an "instantaneous energy loss" . To that end, we take dif- 
ference between the maximum kinetic energy before the 
impact and the next maximum after the impact. Obvi- 
ously the instantaneous energy transfer is practically the 
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same for Nae and Nag independent of their differences in 
structure; and the energy loss is approximately propor- 
tional to E^ in . About half (more precisely around 55%) 
of the impact energy is withdrawn from the cluster in 
that first round. 



3. Redistribution of initial kinetic energy 
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FIG. 14: Distribution of impact energy amongst various com- 
ponents : Kinetic energy of the cluster (split into intrinsic, 
open squares, and translational, close squares, motion), ki- 
netic energy of the MgO substrate (open circles), potential 
energy in the oxygen dipoles (triangles), and other potential 
energy in the substrate (close circles). Upper panel : Nas 
with MgO, lower panel : Na6 with MgO. 

It is, furthermore, interesting to see how the ini- 
tially available energy is distributed over the various con- 
stituents. Such energy balance is shown in figure fT4l as 
a function of E^ in . All energies have been averaged over 
2 ps after impact time. They are drawn relative to the 
maximum kinetic energy of the cluster before impact. 
The energy terms do not necessarily sum up to one be- 
cause the interaction between substrate and cluster is 
omitted, as well as the intrinsic potential energy of the 
cluster. Low energies (< 2.72 eV for Nag and < 2.18 eV 
for Nae) represent the regime of deposition. The largest 
amount of energy is used up here for intrinsic cluster 
motion combined with potential energy of the substrate. 
The latter is responsible for the attachment of the clus- 
ter to the surface. The higher energies represent the 
regime of reflection. The share of energies depends here 
on the cluster geometry. For Na6, the translational mo- 
tion takes the lead, whereas the intrinsic motion shrinks 
to almost zero. That complies with the trajectories in 
figure [5] which show that the cluster is repelled from the 
surface with opposite orientation but only weakly per- 



turbed structure. The reflection of the top ion seems to 
proceed independently from the five-fold ring. The ring 
hits first and departs first while the top ion is reflected 
later thus departing behind the ring. 

For both clusters, kinetic and potential energies of the 
MgO behave similar. For high E® in in the reflection 
regime, potential and kinetic energy are equal. The clus- 
ter quickly transfers some momentum to the substrate 
at impact and then disappears. This leaves the sub- 
strate ions oscillating around their equilibrium positions 
in harmonic motion, associated to equipartition between 
kinetic and potential energy. For high E^ in in the soft 
landing domain, the situation is different. The poten- 
tial energy becomes the dominant contribution because 
the cluster is adsorbed and remains in contact with the 
surface. This distorts the surface and leads to a large po- 
tential energy. Polarization energy is dominating in the 
potential energy. But the isolated contribution from the 
oxygen dipoles, also shown in figure \14\ is comparatively 
small. The polarization within the oxygen ions in fact 
remains small as compared to the polarization caused by 
the displacement of O versus Mg, each one carrying a net 
charge of ±2e. 



IV. CONCLUSIONS 

We have analyzed the dynamics of deposition of small 
Na clusters on an MgO surface taking up a well tested hi- 
crachical QM/MM modeling where the cluster electrons 
are treated quantum mechanically by time-dependent 
local-density approximation and the cluster ions as well 
as the substrate atoms by classical molecular mechanics. 
The dynamical polarizability of the substrate atoms is 
taken into account to describe correctly the strong po- 
larization effects in the cluster-material interaction. The 
results are compared to deposition on Ar surface which 
is much softer than MgO. 

Test cases were Na6, which is a strongly oblate cluster 
with a finite dipole momentum, and the well bound and 
highly symmetrical Nag. The general pattern are simi- 
lar : The clusters are very quickly stopped by the sub- 
strate, they transfer a rather small amount of energy to 
the substrate while acquiring strong internal excitation. 
For larger impact energies, the clusters are reflected from 
the surface. This reflection is, of course, inelastic and 
leaves the clusters departing in highly excited intrinsic 
motion. There are, on the other hand, significant differ- 
ences between the two cases, since the cluster geometry 
has a large influence on the dynamics. The main effects 
come from the strong surface corrugation of MgO (001). 
Nag does not match the surface structure and thus ac- 
quires significant lateral motion in contrast to Nag which 
keeps better on a vertical track. Moreover, Na6 transfers 
more energy to the substrate than Nag. 

In comparison to Ar(001) surface, we find a simi- 
lar energy range for deposition and reflection. The de- 
tails, however, differ dramatically. Deposition on Ar(001) 
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transfers most of the energy to the substrate leaving a 
rather mildly excited cluster on the surface while there is 
very little energy transfer to MgO(OOl) and large intrin- 
sic excitation of the cluster. Reflection from Ar(OOl) is 
achieved at the price of severe surface destruction while 
MgO remains intact at the danger that the highly excited 
departing cluster may fragment later on. 

The detailed energy balance differs in the deposition 
and reflection regime. In case of deposit, most energy 
is going to intrinsic cluster excitation and substrate po- 
larization. In case of reflection, there is, of course, more 
translational energy left for the cluster and the substrate 



develops equipartition of kinetic and potential energy re- 
lated to the remaining small, nearly harmonic, oscilla- 
tions. 
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